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SUMMARY 


The  report  summarizes  the  work  of  the  last  year  under  the  general 
heading  "Effect  of  Local  Material  Imperfections  on  the  Buckling  Behavior  or 
Composite  Structural  Elements".  The  report  contains  two  parts.  The  first 
part  deals  with  delamination  buckling  of  pressure-loaded  long  cylindrical 
shells  and  panels,  and  the  second  with  buckling  of  radially-loaded  annular 
plates  of  laminated  construction  with  or  without  radial  stiffeners. 

In  part  A,  the  geometry  is  virtually  isotropic  and  the  emphasis  is  on 
assessing  the  effect  of  delamination  position  and  size  on  the  critical 
pressure.  The  boundary  conditions  on  the  panels  are  either  simply 
supported  or  clamped  at  circumferential  positions.  Note  that  the 
cylindrical  shells  or  panels  are  extremely  long  and  the  delamination 
extends  along  the  entire  length. 

In  part  B  the  circular  annular  plates  are  subjected  to  uniform  radial 
compression  around  the  circumference.  Laminated  and  isotropic  geometries 
are  used.  For  the  isotropic  geometries  the  effect  of  plate  and  stiffener 
rigidities  on  the  buckling  load  are  fully  assessed.  Moreover,  the  effect 
of  boundary  conditions  is  also  studied.  Finally,  for  the  laminated 
geometry,  a  solution  procedure  is  presented,  and  demonstrated  through  a 
laminated  geometry  that  yields  a  primary  state  which  has  a  uniform  state  of 
stress  (symmetric  and  quasi-isotropic  in  extension). 
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BUCKLING  OF  DELAMINATED,  LONG, 
CYLINDRICAL  PANELS  UNDER  PRESSURE. 


E 


George  J.  Slml. tses+  and  Ziqi  Chen+  + 

School  of  Engineering  Science  and  Mechanics 
Georgia  Institute  of  Technology,  Atlanta,  Georgia 


Abstract 


Delami nation  is  one  of  the  basic  defects  inherent  to  laminar  materials. 
The  investigation  of  the  buckling  characteristics  of  delaminated  cylindri¬ 
cal  shells  or  panels,  when  subjected  to  external  pressure,  is  presented 
herein.  The  geometry  is  such  that  it  covers  a  wide  range  of  length  to 
radius  ratios  as  well  as  panels  of  different  widths.  Results  are  presented 
only  for  very  long  cylinders  and  panels.  The  boundaries  are  either  simply 
supported  or  clamped.  Furthermore,  the  material  is  such  that  it  leads  to 
(quasi)  isotropic  laminates  for  all  sections  involved;  the  overall  as  well 
as  the  ones  separated  by  the  delamination.  Finally,  the  geometry  is  free 
of  initial  geometric  imperfections.  Because  of  the  last  two  assumptions,  a 
primary  membrane  state  exists  and  bi.furcational  buckling  is  possible. 
3uckling  load  are  calculated  for  a  wide  range  of  parameters.  The  width  and 
the  through-the-thickness  position  of  delamination  greatly  affect  the 
bifurcation  load. 


1.  Introduction 


Cylindrical  shells  and  panels  are  widely  used  as  primary  structures  In 
several  applications.  These  are  often  subjected  to  destabilizing  loads. 
Therefore,  buckling  is  an  important  failure  mode  and  it  forms  a  fundamental 
consi derat J on  in  the  design  of  such  systems. 

The  advent  of  fiber  reinforced  composite  materials  has  resulted  in  a 
significant  increase  of  their  use  as  a  construction  material,  because  of 
their  many  advantages,  especially  their  high  potential  weight  and  overall 
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costs  savings.  It  was  deemed  necessary  then  to  investigate  the  buckling 
characteristics  of  laminated  cylindrical  shells.  Initially,  the  studies 
were  confined  to  con f  1  gu rati ons ,  which  are  free  of  defects,  such  as 
del  ami nati ons . 

In  1975,  Tennyson  [1]  presented  a  review  of  all  previous  studies  on 
the  problem.  Most  of  these  studies  employed  the  classical  approach  and 
dealt  with  individual  and  combined  application  of  uniform  axial 
compression,  pressure  and  torsion.  Hirano  [2]  investigated  the  buckling  of 
angle-ply  laminated  circular  cylinders  and  obtained  the  best  lamination 
angles  which  give  the  highest  buckling  stress.  There  are  also  some  papers 
about  the  stability  of  triple-layered  anisotropic  cylindrical  shells  [3“5] 
and  sandwich  panels  [6], 

The  stability  of  geometrically  imperfect,  laminated  cylindrical  shells 
were  investigated  by  Simitses,  Shaw  and  Sheinman  C 7  —  1 0 ]  .  The  governing 
equations  f or  the  nonlinear  analysis  of  imperfect,  stiffened,  laminated, 
circular,  cylindrical  thin  shells,  subjected  to  uniform  axial  compression 
and  torsion,  and  supported  in  various  ways,  were  derived  and  presented. 
Two  types  of  formulations  were  developed,  one  (W,  F  formulations)  based  on 
Sonnell-type  nonlinear  kinematic  relations  and  the  other  (U,  V,  W 
Formulation)  based  on  Sanders-type. 

The  buckling  of  laminated  cylindrical  panels  was  studied  in  recent 
years  also.  Zhang  and  Matthews  [11,12]  considered  panels  under  four  kinds 
or  boundary  conditions  subjected  to  the  combination  of  axial  compression 
and  shear  forces.  Two  coupled,  fourth-order  partial  differential  equations 
were  solved  by  the  use  of  multiple  Fourier  series.  Numerical  computations 
were  performed  for  a  number  of  panels  with  different  lay-ups,  different 
curvatures  and  different  materials. 
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W'nitney  [13]  studied  the  buckling  of  anisotropic  cylindrical  panels 
under  arbitrary  combinations  of  axial  load,  Internal  pressure  and  in-pla-  e 
shear  load.  He  used  Donnell-type  equations  in  conjunction  with  Galerkin’s 
method  to  determine  critical  loads. 

Buckling  and  initial  postbuckllng  behavior  of  symmetrically  laminated, 
thin  cross-ply  cylindrical  panels  under  axial  compression  was  investigated 
by  Hui  [14].  He  obtained  closed  form  solutions  for  the  buckling  loads. 
The  initial  asymmetric  postbuckling  behavior  was  demonstrated  by  computing 
the  postbuckling  coeff icients . 

Composite  structures  often  contain  de 1  ami nat i ons  .  Causes  of 
delamination  are  many  and  include  tool  drops,  bird  strikes,  runway  debris 
hits,  and  manufacturing  defects.  The  presence  of  delamination  in  a 
composite  material  may  cause  local  buckling  and  therefore  a  reduction  in 
the  overall  load  bearing  capacity  of  the  structure.  The  problem  of 
delamination  buckling  has  received  attention  in  recent  years. 

A  finite  element  analysis  was  developed  by  Whitcomb  [15]  to  analyze  a 
laminated  plate  with  a  through-the-wi  dth  delamination.  The  postbuckling 
behavior  was  studied.  In  the  parametric  study,  stress  distributions  and 
strai  n-enery  release  rates  were  calculated  for  various  delamination 
lengths,  delamination  depths,  applied  loads,  and  lateral  deflections.  Some 
delamination  growth  data  were  obtained  through  fatigue  tests.  Another 
paper  in  this  subject  was  presented  by  the  above  author  and  Shivakumar 
[26],  in  1985,  in  which  the  buckling  of  an  elliptic  delamination  embedded 
near  the  surface  of  a  thick  quasi-isotroplc  laminate  was  studied.  Both  the 
Finite  Element  and  the  Rayleigh-Ritz  methods  were  used  for  the  analysis. 
The  Rayleigh-Ritz  method  was  found  to  be  simple,  inexpensive,  and  accurate, 
except  ror  highly  anisotropic  delaminated  regions.  In  that  paper,  effects 
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of  delamination  shape  and  orientation,  material  anisotropy,  and  layup  on 
buckling  strains  were  examined. 

Yin  and  Wang  [17]  derived  a  simple  expression  for  the  energy  release 
rate  associated  with  the  growth  of  a  general  one-dimensional  delamination. 
The  energy  release  rate  was  evaluated  by  means  of  the  path- independent 
J-integral.  Yin  and  Fei  [28]  investigated  the  buckling  load  of  a  circular 
plate  with  a  concentric  delamination. 

Angle-ply  composite  sandwich  beams  with  through-the-width  delami¬ 
nations  were  studied  by  Gillespie  and  Pipes  [19].  Reduction  in  flexural 
strength  was  found  to  be  directly  proportional  to  the  length  of 
delami nation  and  varied  from  41  %  to  87?  of  the  pristine  value. 

S.  S.  Wang  [20]  investigated  the  buckling  of  angle-ply  composite 
laminates  with  edge  delamination.  Based  on  a  recently  developed  theory  of 
laminated  anisotropic  elasticity,  the  problem  was  formulated  using 
'..ekhni  tski  i  '  s  complex  variable  stress  potentials.  An  eigenfunction 
expansion  method  was  employed  to  solve  the  singular  elasticity  problem. 
With  the  aid  of  a  boundary  collocation  technique,  complete  stress  and 
displacement  fields  were  obtained. 

A  two-dimensional  analytical  model  was  developed  by  Chai  and  Babcock 
[21]  t--  assess  the  compressive  strength  of  near-surface  interlaminar 
defects  in  laminated  composites.  The  postbucking  solution  for  the 
delaminated  elliptic  sections  was  obtained  by  using  the  Raylei gh-R i t z 
method,  while  an  energy  balance  criterion  based  on  a  self-similar  disbond 
growth  governed  fracture. 

Simitses,  Sallam  and  Yin  [22-2 iJ  ]  investigated  the  delamination 
buckling  and  growth  of  flat  composite  structural  elements.  A  simple 
one-dimensional  model  was  developed  to  predict  critical  loads  for 
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delaminated  plates  with  both  simply  supported  and  clamped  ends.  The 
effects  of  delamination  poslton,  size,  and  thickness  on  the  critical  loads 
were  studied  in  detail.  The  postbuckling  behavior  as  well  as  the  energy 
release  rate  were  examined.  The  results  revealed  that  the  damage  tolerance 
of  the  laminate  was  either  governed  by  buckling  or  by  the  fracture 
*■  "ugliness  of  the  material. 

Almost  all  the  papers  about  delamination  buckling  deal  with  beams  and 
plates.  Owing  to  its  complexity  in  mathematics,  very  limited  information 
on  the  subject  or  delamination  buckling  of  shells  is  currently  available. 

Troshin  [25]  studied  the  effect  of  longitudinal  delamination,  in  a 
laminar  cylindrical  shell,  on  the  critical  external  pressure.  The  shell 
was  assumed  to  be  separated  by  the  delamination  into  three  panels.  A 
system  or  eight  ordinary  differential  equations  were  derived  from  the 
governing  partial  differential  equations.  The  system  along  with  boundary 
and  continuity  conditions  was  integrated  by  the  Kutta-Merson  method  with 
intermediate  orthonormalization  of  solution  vectors.  Critical  pressures 
ror  various  locations  and  sizes  of  the  delamination  were  found.  In  another 
paper  [26],  the  above  author  Investigated  the  delamination  stability  of 
tr i pie-layered  shells  with  almost  the  same  method.  Sallam  and  Simitses 
[27]  studied  delaminations  buckling  of  thin  cylindrical  shells  of  perfect 
geometry,  when  subjected  to  uniform  axial  compression.  The  delamination 
region  was  assumed  to  be  of  constant  width  and  covering  the  entire 
ci rcumf erence . 

2.  Mathematical  Formulation 

The  Koi  ter-3udi ansky  [28]  buckli  ng  equations  have  been  deduced  from 
those  given  in  the  Appendix  of  [28].  They  are  given  below  in  terms  of 
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stress ,  Njj,  and  moment,  Mjj,  resultants  for  a  thin,  circular,  cylindrical 
shell  under  uniform  pressure. 

M  +  N  -l-M*  +qR(ly°-^)*+qX=o 

xx, x  xy,y  2R  xy,y  4  2  xy  z  ,y 
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where  ;x,  iy  and  dz  are  rotations  about  axes  denoted  by  the  subscripts,  qx, 
q^  and  qz  are  corrections  to  surface  loading  due  to  load  behavior  during 

buckling,  q  is  the  applied  pressure  f°  ,  e°  and  are  the  reference 

xx  yy  xy 

surface  small  additional  strains  (needed  to  take  us  from  the  primary 


"w  *S 
%  * 

fir' 


membrane  state  to  the  infinitesimally  close  buckled  state),  and  R  is  the 
radius  of  the  reference  cylindrical  surface.  Note  also  that  if  the  terms 
denoted  by  a  single  asterisk  are  dropped  the  equations,  Eqs.  (1), 
correspond  to  those  obtained  from  Sanders  type  of  kinematic  relations  [29]. 
Similarly,  if  terms  marked  by  either  single  or  double  asterisk  are  dropped 
one  obtains  the  well  known  Donnell-type  of  equations. 

Since  there  exist  three  possibilities  of  load  behavior  during  the 
buckling  process  [30,31],  the  corrections  to  surface  loading  assume  three 
distinctly  different  expressions. 
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I.  Load  remains  normal  to  the  deformed  surface 
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x  y  ,  v  ,  z 

q  =  -qw,x  ;  qJ  -  ~q  (w,  -  ^ —  )  ;  q  =  0 


(2a) 


II.  Load  remains  constant-directional 


x  y  z  . 
q  =  qJ  =  q  =0 


(2b) 


III.  Load  remains  directed  towards  initial  center  of  curvature 
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The  kinematic  relations  are  given  by 
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where  the  expressions  for  the  reference  surface  strains,  rotations,  <j>j  ,  and 
changes  in  curvature  and  torsion,  kjj,  are  given  by 
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The  usual  lamination  theory  [32]  is  employed.  Moreover,  it  is  assumed 
that  all  shell  (laminates)  parts  are  symmetric  with  respect  to  their 
m j  dsurf aces  (Bjj  =  0).  Then,  the  stress  and  moment  resultants  are  related 
to  the  reference  surface  strain  parameters  by 


k=1 


and  hk  and  hk_-;  denote  the  z-coordinate  of  the  upper  and  lower  surfaces  of 
the  kth  lamina  in  each  laminate,  respectively. 

In  terms  of  the  displacement  components,  u,v  and  w  the  buckling 
equations,  Eqs.  (1),  become 
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Note  that  these  equations  correspond  to  the  Sanders'  approximation  [30]  and 
therefore  all  asterisks  have  been  dropped. 

2 . i  Description  of  the  Delamination  Problem 

Thin  circular  cylindrical  shells  and  panels  with  longitudinal 
delamination  over  the  entire  length  are  next  considered.  The  geometry, 
loading  and  coordinate  systems  are  shown  on  Fig.  1.  The  ends  of  the  panel 
are  either  clamped  or  simply  supported.  The  location  and  size  of  the 
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del  anl  nation  i.s  arbitrary.  Angle  a  denotes  the  region  of  the  delamination , 
while  B  and  v  denote  the  location  of  it  from  left  end  and  right  end, 
res pec ti vely . 

It  is  assumed  that  under  subcritical  loads  the  delamination  does  not 
expand.  The  panel  is  separated  into  four  parts  (four  panels)  by  the 
del  ami  nation.  Each  part  has  a  set  of  coordinates  attached  to  it,  (see  Fig. 
1.)  and  the  natural  plane  of  the  panel  lies  on  the  XY-plane.  The  panel 
is  subjected  to  uniform  external  pressure,  q,  over  the  entire  outer  surface. 
Let  h’  (i  =  I,  II,  III,  IV)  denote  the  thickness  of  the  ith  panel  (see  Fig. 

l).  The  nondimensional  parameter  h  =  h^/h^1'  is  used  to  describe  the 
thickness  of  the  delaminati  on.  Let  u* ,  v* ,  wJ  (i  =  I,  II,  III,  IV)  be  the 
displacement  components  of  material  points  on  the  midplane  of  each  part 
(each  panel)  in  the  x,  y,  and  z  directions,  respectively. 

The  panel  becomes  a  complete  cylindrical  shell  when  the  total  angle 
:(:  =  c  +  b  +  y)  equals  2v .  The  geometry  is  such  that  a  membrane  primary 

state  exists  (3^  3  0)  for  all  participating  parts.  Therefore,  the  buck¬ 
ling  equations,  Eqs.  (1),  apply  to  each  part.  They  are  subject  to  boundary 
conditions  at  03  =  0  and  64  =  y  (see  Fig.  1),  and  certain  kinematic  con¬ 
tinuity  conditions  as  well  as  force  and  moment  local  equilibrium  conditions 
at  the  common  boundaries  of  the  various  parts  (63  =  B  and  04  =  0). 

i.'ote  that  all  of  these  conditions  are  associated  with  y  =  constant 
positions.  Each  part  must  also  satisfy  boundary  conditions  associated  with 
x  =  constant  positions.  These  are  not  listed  now,  but  they  correspond  to 
the  classical  simply  supported  ones.  They  are  listed  in  the  next  section. 

The  various  boundary  and  auxiliary  conditions  associated  with  y  = 
constant  positions,  are  listed  below: 
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Auxiliary  Condition  at  83  =  8  (9i  a  8; 
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Auxiliary  Conditons  at  8a  =  0  (61  =  6?  =  a) 
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Boundary  Condition?  at  9t  =  I 
(a)  Clapped 

wIV  =  0 


(b)  Simply  Supported 
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Note  that  the  simply  supported  conditions,  Eqs.  (8b)  and  (11b) 
correspond  to  the  classical  simply  supported  conditions,  SS~3  [32]. 


L2 


3.  Solution  Procedure 
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where  the  prime  denotes  differentiation  with  respect  to  0, 


L  =  L/rrmR  ,  0  =  y/R,  and  A  =  qR  /D  , 


and  R  is  the  same  for  all  parts. 

Elimination  or  U  and  V  through  the  use  of  the  first  two  of  Eqs.  (15) 
and  through  substitution  into  the  third  one  yields  a  single  higher  order 


!  j nary 


differential  equation  in  W  alone.  This  higher  order  equation 


assumes  the  form 


. • .  - .. 


where  the  F's  are  constants  that  contain  the  external  load  q  and  structural 


geometri  c  parameters  (Ajj,  Djj,h,  R,  etc.).  Note  that  some  of  the  F's 
change  according  to  the  case  of  load  behavior  during  buckling,  Eqs.  (2), 
The  expressions  for  the  F's  are  given  below: 
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e  colu.T.n  vector  in  the  expression?  for  ajj  denote  the  load  cases  as 

ai  h>  F  :s .  '2). 


The  solution  procedure  Is  similar  to  the  one  described  for  rings  and 
arches  in  Article  7.3  of  Ref.  31).  The  number  of  equations  is  higher,  the 
equations  themselves  are  more  complex,  and  a  closed  formed  solution  is  not 
expected  as  in  Ref.  33.  Nevertheless,  the  overall  procedure  can  be 
followed  and  a  numerical  estimate  can  be  achieved. 

The  basic  steps  are  as  follows: 

Fjrst  assume  for  W(9)  a  solution  of  the  form 


W  =  C  er9  (20) 

Since  the  order  of  the  equation  is  eight,  then  substitution  into  Eq. 
(17)  yields  an  eighth  degree  polynomial  in  r.  Thus,  eight  roots  are 
expected  for  each  geometry  and  load  level.  If  the  eight  roots  are 
distinctly  different,  the  general  solution  for  W(9)  is  given  by 

8  r .  6 

W(9)  =  I  C  e  3  (21 ) 

i=1 

Note  that  if  double  roots  occur  the  form  of  Eq.  (21)  is  modified 
accord! ngly . 

The  form  of  the  solution  for  U(B)  and  V(e)  is  similar  to  that  of  Eq. 
( 21  ) ,  or 
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(22) 


Note  that  in  deriving  Eq.  (17),  U  and  V  were  eliminated  through 
operations  on  Eos.  (15).  The  intermediate  steps  that  lead  to  this 


elimination  yield 


(23) 


(L21L12  ”  L11L22)  U  (L13L22  “  L 1 2L23 ^  W 


(L21L12  "  L1 1L22)  V  (L11L23  ~  L13L21*  W 


Subs  t  j  t  ut  j.  on  of  Eqs.  (21)  and  (22)  into  Eqs.  (23)  yields  the 
expressions  of  and  Bj  in  terms  of  Cj,  or 

A  -  C-C. 
j  s  j  j 


B.  =  n.C. 
).  J.  J. 


for  i  =  1 , 2 , 8 


where 


(a11  +  ria12  +  ria13)/(a31  +  ria32  4  ria33) 


(a2iri  +  ria22  4  ria23)/(a31  +  ri332  +  ria33) 


(24a) 


Tnus,  the  solutic.  for  U,  V,  and  W  is  given  in  terms  of  eight  constants  Cj, 
while  and  n*  are  known. 


U  (9)  =  l  $  C  e 
i=1 


3  r  .9 

V  (9)  =  I  n . C j  e 


8  r .  0 
W  (0)  =  J  C  e  5 


Tnere  exist  eight  unknowns,  Cj's,  for  each  panel.  Since  there  exist 
four  panels  be  'ause  of  the  delamination,  the  total  number  of  unknowns  is  32. 
ire  or  the  32  boundary  and  auxiliary  conditions,  Eqs.  (8)— (11),  leads  to  a 
system  32  linear  homogeneous  algebraic  equations  in  the  32  unknowns.  A 


nontrivial  solution  exists  if  the  determinant  of  the  coefficients  vanishes. 


This  i s  the  characteristic  equation,  which  yields  the  value  of  the  critical 
load.  A  numerical  procedure  has  been  developed  to  accomplish  this. 


For  the  sake  of  completeness,  the  boundary  and  auxiliary  conditions 
are  presented  in  terms  of  the  respective  displacement  functions,  u,  v  and  w 
[see  Fqs.  ( 8 ) - ( 1 1 )]. 
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Boundary  Conditions  at  9a  -  7 


(a)  Clamped 


u  =  0 


( b )  Simply  Suppnr ted 
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All  of  the  above  conditions  can  easily  be  expressed  in  terms  of  U(e), 
7(9),  and  W(9)  instead  of  u,  v,  and  w.  This  is  so  because  in  all  boundary 
and  auxiliary  conditions,  Eqs.  (26)-(29),  if  the  order  of  differentiation 
with  respect  to  x,  of  v  and  w  is  even,  the  one  for  u  is  odd,  or  vice  versa 

[see  Eq.  (19)],  Thus,  in  each  one  of  the  conditions  either  si  n  or 

Li 

'T’lmy  x 

cos  is  the  common  factor,  which  does  not  vanish  for  all  x. 

A  computer  program  has  been  written  in  order  to  obtain  critical 
conditions  f or  all  geometries.  The  Georgia  Tech  high  speed  digital 


computer,  CDC  Cyber  70,  Model  79-28,  was  used  for  generation  of  results. 


4.  Numerical  Result? 


Results  are  presented  herein  only  for  a  special  geometry  and  only  load 
case  I,  Eqs.  (2a). 

The  geometry  Is  an  Isotropic  geometry  (as  In  [25])  w  j  t  h  the 
delami nation  parallel  to  the  shell  reference  surface  and  extending  along 
the  entire  length  of  the  panel.  Moreover,  the  panel  is  assumed  to  be  very 
long  and  thus  the  solution  is  not  affected  by  the  x-constant  boundary 
conditions.  Finally,  the  boundaries  at  6  =  const,  are  taken  to  be  clamped. 
The  results  are  used  to  study  the  effect  of  delamination  position 

(through  the  thickness),  h,  and  of  delamination  size,  a  (see  Fig.  1).  Note 
that  the  results  correspond  to  0  =  Y  for  the  case  of  panels. 

The  results  are  presented  both  in  tabular  and  graphical  form,  for  a 
complete  circular  shell  as  well  as  for  panels  of  angle  ir  and  tt/2.  Critical 
load  parameter  values,  -A  (1  =  qR3/D),  are  shown  in  Tables  1~3  for  various 

values  cr  3  and  h.  The  mode  is  either  symmetric  (S)  or  antisymmetric  (A), 
and  it  is  designated  as  such  in  the  tables  and  on  the  figures,  Figs.  2-4. 

For  the  complete  cylindrical  shell,  and  for  midsurface  delamination 

n  =  0.5,  the  critical  load  parameter  is  3  for  the  perfect  configuration, 
and  it  decreases  to  0.75  when  the  delamination  extends  to  the  entire 
ci  rcu.m^erence .  It  is  next  shown  that  this  is  a  reasonable  expectation. 

It  is  assumed  that  before  buckling  the  circularity  is  maintained, 
there  is  complete  contact  of  parts  I  and  II  (see  Fig.  1),  and  if  one 
ienctes  by  r:} ,  and  c}  ‘  the  pressures  on  the  various  parts  then 


is 


an  : 


(30) 


Moreover,  in  the  primary  state 
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and  R3  =  R 


Then,  before  and  up  to  the  instant  of  buckling  the  pressure  loadings  on 
parts  I  and  II,  q1  and  q11,  are  related  to  the  applied  loading  q,  through 
Eqs  .  ( 31 ) ,  by 


q  -  h  q 


q  =  (1  -  h)  q 


Ir  A  is  defined  as 

A1  =  q5R3/  D1 


1  =  I,  II,  III,  and  IV 
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Clearly  then,  when  the  cylindrical  shell  is  completely  delaminated 

A  =  h2AX 
cr  cr 

and  A1  =  ~3.0;  therefore  A  =  -  0.75  for  h  =  0.5. 
cr  cr 

One  may  say  at  this  point  that  the  residual  strength  of  the  completely 
delaminated  thin  cylindrical  shell  is  equal  to 


Thus,  for  h  =  0.5 


1  total  '  '  C-S)2](-3) 

A,  .  ,  -  “1.5,  while  for  h  =  0.1 
total 


' total  =  0.82(-3)  =  “2.46. 


The  results  corresponding  to  the  complete  thin  delaminated  cylinder 
are  shown  in  Table  1.  Note  that  when  a  =  6.2  (close  to  2tt)  the  value  of 
the  buckling  load  parameter  is  (0.1)2  (3)  or  0.0003.  The  same  results  are 
presented  graphically  on  Fig.  2,  while  Fig.  3  shows  the  complete  curves 

corresponding  to  symmetric  and  antisymmetric  buckling  but  only  for  h  =  0.5. 
Furthermore,  there  exists  two  important  observations.  One  is  that  the 

critical  load  associated  with  the  very  small  h-values  should  be  viewed  as 
measures  0^  local  buckling  and  not  as  measures  of  load  bearing  capacity  of 
the  overall  structure  [25].  The  second  observation  is  that  in  most  cases 
the  buckling  mode  of  the  two  parts  that  are  separated  by  the  delamination 
(parts  I  and  II  ;  see  Fig.  1),  is  an  antisymmetric  mode  of  the  type  that 
suggests  contact  over  a  certain  portion  and  no  contact  over  the  remaining 

one.  This  suggests  that  the  postbuckllng  behavior  for  shells  with  h  will 
be  di  Cerent  from  that  of  shells  with  (1  -  h),  although  the  critical  loads 
are  the  same  for  the  two  (see  Table  1).  Moreover,  for  each  geometry, 

regardless  of  h,  the  transition  from  buckling  to  postbuckllng  behavior  will 
require  accommodation  of  the  predicted  (buckling  mode)  contact. 

Results  for  clamped  panels  with  <i>  =  tt  and  tt/2  and  symmetric 
delamination  (6  =  "v )  are  presented  in  Tables  2  and  3  and  on  Figs.  H  and  5. 
All  generated  results  correspond  to  an  antisymmetric  mode.  For  these 
results  also,  as  in  the  case  of  the  complete  cylinder,  the  critical  load 
varies  rrom  the  perfect  geometry  critical  load  [3*1,35],  when  a  is  very 

-2  -  2 

small,  to  the  value  of  h  Aperf.  [°r  (i-1"1)  A  perf.],  when  the  panel  is 
completely  delaminated  (see  Fig.  1). 


From  the  above  results  it  is  clear  that  further  and  more  detailed 


studies  need  be  performed  before  we  fully  understand  the  complex  response 
op  delaminated,  curved  laminates. 
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Table  1  Buckling  Load  for  a  Delaminated  Circular  Cylindrical  Shell 


0.5 

0.3(0. 7) 

0.1(0. 9) 

0.01(0.99) 

0.000 

3.00000(A) 

3.00000(A) 

3.00000(A) 

3.00000(A) 

0.2 

2.88461  (S) 

2.86341 (S) 

2.93602(S) 

0.20186 (A) 

0.4 

2.41778 (S) 

2.57496 (S) 

2. 87696 (S) 

0.05044(A) 

0.6 

2.27045  (S) 

2. 46545 (S) 

2.23770(A) 

0.02240(A) 

0.8 

2.23305  (S) 

2.42099 (S) 

1.25739 (A) 

0.01260(A) 

1.0 

2.19743(A) 

2.40772 (S) 

0.80364(A) 

0.00803(A) 

1.2 

1.94698(A) 

2.31670(A) 

0.55715(A) 

0.00557(A) 

1.4 

1.76893(A) 

2.12619(A) 

0.40858(A) 

0.00408 (A) 

1.6 

1.65424(A) 

1.95652(A) 

0.31218(A) 

0.00312(A) 

1.8 

1.58796 (A) 

1.78177(A) 

0.24614(A) 

0.00246 (A) 

2.0 

1.55616(A) 

1.57452(A) 

0.19894(A) 

0.00199(A) 

2.2 

1.54597(A) 

1.35834(A) 

0.16406 (A) 

0.00164(A) 

2.4 

1. 50758 (S) 

1.16633(A) 

0.13757(A) 

0.00137(A) 

2.6 

1. 38708 (S) 

1.00618(A) 

0.11701(A) 

0.00117(A) 

2.8 

1 . 291 18 (S) 

0.87462(A) 

0.10075(A) 

0.00101(A) 

3.0 

1 .21859 (S) 

0.76648 (A) 

0.08769(A) 

0.00088(A) 

3.2 

1. 16706 (S) 

0.67713(A) 

0.07706(A) 

0.00077(A) 

3.4 

1.13386(5) 

0.60284(A) 

0.06832(A) 

0.00068 (A) 

3.6 

1.11579 (S) 

0.54070(A) 

0.06107(A) 

0.00061(A) 

3.8 

1 . 10883 (S) 

0.48847(A) 

0.05501(A) 

0.00055(A) 

4.0 

1. 10645 (S) 

0.44445(A) 

0.04993(A) 

0.00050(A) 

4.2 

1.03733(A) 

0.40729(A) 

0.04566(A) 

0.000 46(A) 

4.4 

:  0.97570(A) 

0.37596(A) 

0.04207(A) 

0.00042(A) 

4.6 

0.92187(A) 

0.34970(A) 

0.03907(A) 

0.00039(A) 

4.8 

0.87585(A) 

0.32788(A) 

0.03658(A) 

0.00037(A) 

5.0 

0.83751(A) 

0.31006 (A) 

0.03455(A) 

0.00035(A) 

5.2 

0.80665(A) 

0.29590(A) 

0.03294(A) 

0.00033(A) 

5.4 

0.78308(A) 

0.28514(A) 

0.03172(A) 

0.00032(A) 

5.6 

0.76645(A) 

0.27756(A) 

0.03086(A) 

0.00031(A) 

5.8 

0.75620(A) 

0.27287(A) 

0.03032(A) 

0.00030(A) 

6.0 

0.75131(A) 

0.27061 (A) 

0.03007(A) 

0.00030(A) 

6.2 

0.75003(A) 

0.27002(A) 

0.03000(A) 

0.00030(A) 

SYMM.  Buckling  Mode 
ANTI-SYMM.  Buckling  Mode 


Table  2  Buckling  Load  for  a  Delaminated  Panel  of  Angle  7f(rad.) 
_ (Clamped  Boundaries) 


0.5 

0.3(0. 7) 

0. 1(0.9) 

0.01 (0.99) 

0.000 

8.00000(A) 

8.00000(A) 

8.00000(A) 

8.00000(A) 

0.2 

7.90294(A) 

7.94443(A) 

7.98782(A) 

0.20187(A) 

0.4 

7.18972(A) 

7.56860(A) 

5.03907(A) 

0.05044(A) 

0.6 

5.81200(A) 

6.71586(A) 

2.23788(A) 

0.02240(A) 

0.8 

4.64551 (A) 

5.71862(A) 

1.25744(A) 

0.01258(A) 

1.0 

3.94345(A) 

4.86852(A) 

0.80366(A) 

0.00804(A) 

1.2 

3.58022(A) 

4.10046 (A) 

0.55720(A) 

0.00557(A) 

1.4 

3.43175(A) 

3.31777(A) 

0.40862(A) 

0.00409 (A) 

1.6 

3.40141(A) 

2.63529 (A) 

0.31223(A) 

0.00312(A) 

1.8 

3.39324(A) 

2.12076(A) 

0.24619(A) 

0.00246 (A) 

2.0 

3.31556(A) 

1.73589 (A) 

0.19898(A) 

0.00199 (A) 

2.2 

3. 13397(A) 

1.44421 (A) 

0.16410(A) 

0.00164(A) 

2.4 

2.88519(A) 

1.21911(A) 

0.13762(A) 

0.00138(A) 

2.6 

2.61981(A) 

1.04247(A) 

0.11706 (A) 

0.00117(A) 

2.8 

2.36815(A) 

0.90161 (A) 

0.10080(A) 

0.00101(A) 

3.0 

2.14240(A) 

0.78775(A) 

0.08774(A) 

0.00088(A) 

3.141 

2.00056(A) 

0.72026(A) 

0.08003(A) 

0.00080(A) 

(S)  SYMM.  Buckling  Mode 

(A)  ANTI-SYMM.  Buckling  Mode 


Table  3  Buckling  Load  for  a  Delaminated  Panel  of  angle  7T/2(rad.) 

(Clamped  Boundaries) 


0.5 

0.3(0. 7) 

0.1(0. 9) 

0.01 (0.99) 

0.000 

32.5951A(a) 

32 . A9695 (A) 

32. A337A (A) 

32. A3112(A) 

0.2 

29.21156(A) 

30.70972(A) 

20.16636 (A) 

0.20186(A) 

0.  A 

18.67281(A) 

23.12066(A) 

5.03951(A) 

0. 050AA (A) 

0.6 

13.92583(A) 

16.33205(A) 

2.23812(A) 

0. 022A0 (A) 

0.8 

12.87370(A) 

10.5A976 (A) 

1.25768 (A) 

0.01258(A) 

1.0 

12.7255A(a) 

6.99228(A) 

0.80391(A) 

0.00803(A) 

1.2 

1 1 . A62A9 (A) 

A. 93227 (A) 

0.557A5 (A) 

0.00558(A) 

1.  A 

9.57106 (A) 

3.65625(A) 

0. A0889 (A) 

0.00A09 (A) 

1.57 

8 . 1 139A (A) 

2.92168(A) 

0. 32A6A (a) 

0.00325(A) 

SYMM.  Buckling  Mode 
ANTI-SYMM.  Buckling  Mode 
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EFFECT  OF  BOUNDARY  CONDITIONS  AND  RIGIDITIES 
ON  THE  BUCKLING  OF  ANNULAR  PLATES 

G.  J.  Simitses  and  Y.  Frostig 

Engineering  Science  and  Mechanics 
Georgia  Institute  of  Technology,  Atlanta,  GA.  30332 

ABSTRACT 

The  buckling  analysis  of  multi-annular  plates,  with  different 
properties,  is  presented.  The  geometry  may  include  ring  stiffeners  at  the 
common  joints  and  the  plate  may  be  supported  in  various  manners,  simple 
supports,  clamped  supports,  etc.,  at  the  loaded  outer  edge.  The  loading  is 
uniform  radial  compression  and  static.  Several  parametric  studies  are 
performed  in  order  to  assess  the  effect  of  geometry,  material  properties  of 
the  annular  sections  and  of  the  ring  stiffness.  Moreover,  when  rings  are 
present,  the  ring  geometry  is  modelled  both  as  a  curved  beam  and  as  an 
annular  plate,  for  comparison  purposes. 

INTRODUCTION 


The  stability  of  multi-annular  plates  with  or  without  ring  stiffeners 
is  being  considered.  Most  researchers  have  investigated  the  buckling  of  a 
single  annular  plate  (with  or  without  a  rigid  inclusion)  when  subjected  to, 
primarily,  a  uniform  stress  field. 

Several  studies  [1-8]  nave  been  reported  in  the  open  literature  for 
circular  and/or  annular  plates  under  various  loads,  boundary  conditions  and 
thickness  variation;  see  also  their  cited  references. 

Moreover,  a  few  studies  report  on  stiffened  configurations  with 
special  stiffening  [4,9,10]. 

In  all  of  the  above  studies,  certain  simplifying  assumptions  were 
made,  such  as  neglecting  the  extensional  stiffness  of  the  ring  stiffeners 
[10],  which  nay  or  may  not  have  a  significant  effect  on  the  critical  load. 
Moreover,  the  annular  sections  were  considered  to  be  homogeneous  and  of  the 
same  material,  by  virtually  all  investigators. 

The  present  manuscript  deals  with  buckling  of  multi-annular  plates 
with  different  material  properties  and  thickness  and  stiffened  or 
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ans  t  i  f  r'ened  by  rings  at  the  common  joints.  The  inplane  loading  is 
a xi symmetric  and  several  combinations  of  transverse  and  inplane  boundary 


conditions  are  considered.  The  mathematical  modeling  accounts  for  the 


extensional  rigidity  of  the  rings,  in  addition  to  their  torsional  and 


bending  stiffnesses. 


MATHEMATICAL  FORMULATION 


The  equations  governing  the  primary  state  for  each  part  are;  (see 


■'r'J  )  -  N 

'  rr;’r  60 


N  =  0 

00,  6 


Tr.e  buckling  equation  for  each  part  of  the  plate  is: 

(2> 

r  r 

0  o 

where  N  ,  are  primary  state  stress  resultants  in  the  radial  and  cir- 
rr  ee 

oumferential  directions,  respectively;  w  is  the  transverse  displacement; 

(  ),  (  ),  are  partial  derivatives  with  respect  to  radial  or  angular 

P  O 

coordinates  and  D  is  the  flexural  rigidity.  The  geometry  and  loading  are 
sncwn  on  figure  1. 

The  inplane  stress  resultant  distribution  is  solved  first,  using 
equations  (1).  The  constants  are  determined  through  continuity  conditions 
at  tne  common  joints,  see  [11].  In  case  of  a  stiffened  plate,  the 
continuity  conditions  include  the  deformation  and  radial  stress  resultants 
f  the  stiffener  (see  [12]). 

Tr.e  ir.ulane  stress  distribution  in  the  various  parts  of  the  plate 
ce; ends  on  the  extensional  rigidity  of  the  various  parts,  in  addition  to 
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the  extensional  rigidity  of  the  stiffener.  Even  in  case:;  where  r  :  :  .  > ... 

is  homogeneous,  but  with  stiffeners,  the  stress  distributi  .■  r:  is  r:d 
constant  as  suggested  by  Rossettos  [10].  The  exact  distribution  defends  n 
the  extensional  rigidity  and  location  of  stiffeners.  The  buckling  equation 
for  each  part  is  given  by  equation  (2).  The  equation  is  solved  by  assuring 
a  separable  solution  of  the  form: 

so 

w  ( r  ,  0 )  =  I  F(r)  cos  no  ;  o'.- 

n  n 
n=0 

wnere  Fn(r)  is  the  deflection  function  in  the  r  direction  for  the  r.th 
buckl ing  node . 

The  constants  for  the  transverse  deflection  are  calculated  using 
continuity  and  equilibrium  conditions  at  the  common  joints  and  boundary 
conditions  at  the  edges;  see  [11]  for  the  unstiffened  plate.  In  the  case 
of  a  stiffened  plate  the  continuity  conditions  are  applied  also  to  the 
stiffener.  The  equilibrium  conditions  include  the  contr  *  out  i  or.s  of  the 
stiffener,  related  to  their  bending  and  torsional  rigidities;  sue  [12]. 

In  the  case  of  a  plate  with  a  ring  of  rectangular  cross-sect i on  ,  it  is 
possible  to  analy-e  it  by  considering  the  entire  system  as  a 
mul ti-annnuilar  plate.  In  such  a  case  the  dimensions  of  the  rectangular 
ring  must  be  such  that  plate  theory  is  applicable  to  all  plate  segments 
(inner,  ring  and  outer).  Here,  the  various  ring  stiffnesses  car.  be 
expressed  in  terms  of  the  ring  plate  material  properties  ( K;..  ,  v 
thickness  tp .  in  this  case,  if  the  bending  rigidity  i  a  f  i  «  i  ,  l  r 
torsional  rigidity  is  also  fixed  (do;  ending  only  or.  the  Psis:.  s  i  o:  ;  . ) . 

A  numerical  example  is  presented  i  r.  a  sutse  puent  seel  in  t 

results  are  compared  witr.  these  based  or.  r- i  ng-st  i  ft  or.e  :  ;  1  .it*  or.;'...  . .. . 

-K" 


s , 


% 
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NUMERICAL  EXAMPLES 


A  computer  program  has  been  developed,  and  it  is  described  with 
sufficient  detail  in  [11].  With  the  aid  of  this  program  several  parametric 
studies  are  performed.  The  results  consist  of  critical  loads, 
corresponding  to  axisyrnmetric  (n  =  0)  and  asymmetric  (n  =  1 )  buckling 
modes,  for  multi-annular  plates.  The  geometry  is  a  two  part  plate,  one 
annular  and  one  circular,  with  different  material  properteis,  with  or 
without  a  single  ring  stiffener,  located  either  at  the  edges  or  at  the 
common  joint. 

The  stiffener  rigidities  are  expressed  in  a  non-dimensional  form  as 
follows: 

=  SAft  .  3  =  _ _ 

’  etV1  -  ^  ro  ’  R  Vi712  (1  -  ro 

Yr  =  CJR/  EIR  W 

wr.ere,  EA^>  EIr,  GJrt  are  the  extensional,  flexural  and  torsional  rigidities 
of  the  ring,  respectively.  E1,t1,v1  are  the  modulus  of  elasticity, 
thickness  and  Poisson's  ratio,  respectively,  for  the  outer  plate,  and  r0  is 
tr.e  outer  radius  of  plate. 

3efore  discussing  the  numerical  results,  we  present  the  influence  of 
the  extensional  rigidites  of  the  various  sections,  including  the  one  of  the 
stiffener,  on  the  distribution  of  the  primary  (prebuckling)  stress 

..  ...  ... 


Primary  State  Distributions 


The  primary  state  is  solved  for  first,  by  equations  (1).  The  presence 
of  a  stiffener,  because  of  its  extensional  rigidity,  influences  the  stress 
di str  i  but  i  ons .  The  parameters  to  be  considered  are  the  ratio  of  the 
moduli,  E2/^i,  the  stiffener  rigidity,  ctR,  and  the  position  of  common 
joint,  r i /r0 . 

The  constant  stress  resultant  of  the  inner  plate  is  derived 
analytically,  by  using  equations  (1)  and  the  continuity  conditions. 


)  1 


N° 
rr , 


J0  J  l.f-  t)(.  .  v) 


rr 


(1 


v)  B1  _ 


2  l 


)  ] 


(5) 


where  u./u  =  2 
1  o 


B 1  /  [  (1  ♦  v)  +  B*  O  -  v)  «•  (E2/E1)  (1  +  v)  (B*  ~  D 


aR61  (Bf  ~  D  ] 


3i  =  ro/rr 


(6) 


v  -  Poisson's  ratio  of  the  plates  (v1  *>  v2  =  v) 


and  N 

rr 


applied  at  the  outer  edge. 

In  case  of  a  plate  stiffened  by  an  edge  ring,  the  constant  stress 
resultant  of  the  inner  part  becomes 


N°  =  -  M  /  [  1  ♦  cuE  /E  (1  +  v)]  (7) 

rr  ^  rr  R  1  2 

The  effect  of  the  stiffener  extensional  rigidity  on  the  stress 


resultant  distribution  of  an  annular  plate,  stiffened  by  a  single  ring  at 
tne  inner  edge,  appears  on  figure  2.  As  ap  is  increased  the  stress 


(5  throughout  the  annular  plate.  Note  that  if  there  is  no 
r.e  radial  stress  resultant  is  zero  and  increases 
o  the  value  of  the  applied  loading,  Nrr.  When  ap  *  0. 

once  of  the  ring,  the  value  of  N°  close  to  the 

rr 

ter  than  zero  _nd  its  variation  with  position  (r/rG)  is 

if  a  multi-annular  plate,  stiffened  by  a  ring,  the  stress 
ion  depends  on  ap  and  ring  position  r^/r0,  but  also  cn 
f  tr.e  two  plate  sections,  E2/E-|  -  Figure  3  describes  the 
for  a  two-part  plate  stiffened  by  a  ring  at  the  common 
|/r0  =  0.7).  The  ring  extensional  stiffness  parameter  is 
( corresponding  to  light  stiffening  L12])  and  E2/E-|  = 

.  Note  that  for  E2/Ej  =  1.0,  the  stress  resultant 
e  towards  the  ring,  then  there  is  a  drop  (because  of  the 
ess  cf  the  ring),  and  finally  it  becomes  constant  for  the 
'2/1;  =  0.5  and  0.1,  the  only  difference  is  that  the 
listribution  cecreases  (outer  annular  section)  as  we  move 

It i- annular  Plates 

plate  is  s i mpl y- s uppor t ed  and  loaded  by  a  radial 
ie  loading  at  the  outer  edge.  The  ease  of  a  stiffened 
n  a  single  ring  located  at  the  common  joint,  is  also 

are  calculated  for  the  following  parameters: 

;  n  -  0,  1  (modes);  ti  ■=  t2; 

’ l o ) ,  0.01,  0.1,  0.5,  1.0; 

.0,  0.1  .  0.9,  1.0. 
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<:-r.eu  c  >nf  iguratior.  the  following  parameters  are  considered: 

=  1.2  (axisym.) 

X\>,  :  3-dc;  Sr  =  1.2;  Yr  =0.6  (asym.) 
torsional  ring  stiffness,  Yr,  does  not  affect  axisymmetric 
ing. 

alts  for  the  unstiffened  plate  appear  on  figure  4  for  the 
buckling  mode.  For  the  same  mode,  the  results  for  the 
.-'.figuration  appear  on  figure  5.  Figures  6  and  7  describe  the 
ospor.ding  to  the  asymmetric  mode  (n  =  1),  for  unstiffened  and 
tea,  respectively. 

owing  conclusions  can  be  drawn  from  these  figures: 

critical  load  (lowest  eigenvalue)  corresponds  to  the 
'.ode,  for  both  types  of  plates. 

buckling  load  for  an  unstiffened  plate  is  decreasing  as  the 
e  radii  increases.,  This  is  true  for  both  modes.  In  the 
3,  figure  5,  the  critical  load  increases  with  r-|/r0,  for  r^/r0 
-.e.n  it  decreases,  for  r-\/rQ  £  0.7.  For  a  stiffened  annular 
-  j,  the  critical  load  increases  continuously  with  increasing 

a  two-part  unstiffened  plate  (0  S  =  1.0),  the  critical 

;n  an_.  f\/r0.  For  extremely  small  values  of  E2/E^  (it 

r  values  of  0.1  and  0.01,  and  the  observation  probably  is  true 
■  va.ue 3  in  the  range  of  0.1  <  E  2  /  E ;  <  0.5  the  behavior 
■  a  ratio  of  the  radii  (r,/r0).  For  E-./E,  =  0.1,  as  long  as 
er  t .-.an  approximately  0.55,  buckling  seems  to  to  triggered 
m-.  annular  part  and  the  critical  load  is  a  little  larger  than 
..■'.Hr.  •/  to  E.?/Ei  =  1  (pure  annular  plate).  ..  n  the  ether 

4-  X 


1  loads  for  a  ring-stiffened  plate 


ri/rQ  -  values  larger  than  0.55,  the  critical  load  is  smaller 
of  the  pure  annular  plate,  and  buckling  is  triggered  by  the  inner 
flexible  part.  It  appears  that  for  E2 /E i  =  0.1,  the  corresponding 
caches  the  value  of  0.4282  as  r^/r 0  approaches  unity  (a  complete 
plate  with  bending  rigidity  equal  to  one  tenth  of  that  of  the  case 

Ir.  the  stiffened  case,  the  central  part  triggers  overall  buckling 
tiffer.er  is  located  between  0.7  <  r  1  / r 0  £  1.0.  For  values  of 
.7,  the  modulus  of  elasticity  ratio,  E2/E1  has  small  influence, 
itical  loads  are  nearly  the  same  as  those  corresponding  to  the 
annular  plate.  For  r^/rc  £  0.7  (stiffened  plate)  the  decrease  in 
eui.  load  depends  on  L'^/Ei  .  As  the  inner  part  becomes  more 
tr.o  "drop"  in  the  critical  load  becomes  much  steeper. 

Fir  2  stiffened  plate  with  an  outer  edge  ring  (see  figure  5),  the 
sad  coes  not  converge  to  the  expected  value  of  4.20  D^/r  ,  but 
or  value.  This  is  so  because  of  the  extensional  rigidity  of  the 
equation  (7)). 

For  0  <  Zp/Ej  £  1.0,  as  the  stiffener  approaches  the  plate  edge, 
a i ring  effect  of  the  ring  is  more  pronounced  and  buckling  occurs 
plate  (inner  part)  is  partially  clamped. 

r,  one  stiffener  geometry  (cir  =  0.08,  Br  =  1.2  and  Yr  =  0.6)  is 
tr.e  study  was  concentrated  on  the  effect  of  properties  of  the  two 
s  C  =- 2 / — 1  l1  and  the  stiffener  position,  r]/rQ,  for  both  modes  (n  = 


the  effect  of  stiffener  rigidities  is  studied  for  two  Eg/E} 

.  ar:i  1.1)  for  both  axisymmetric  and  asymmetric  modes.  It  is 


?,  though,  for  only  one  value  of  or  (  =  0.10).  The  geometries  are 

rr.ed  by: 

n  =  0;  =  0.10;  v  =  0.3;  E2/Ei  =  1.0,  0.1;  and  Br  =  1,2,5,10 

,  axi symmetric) 

r.  =  1  ;  1R  =  0.  10;  v  =  0.3;  E2/E1  =  1.0,  0.1  ;  Br  =  5;  and  Yr  =  0.3, 

:.7,  1.0,  1.5,  2.0. 

Toe  results  are  shown  on  figures  8  and  9  for  n  =  0  and  n  =  1, 
e  at  i  ve  1  y . 

Tr.e  general  trends  are  similar  to  those  described  previously, 
over,  the  ring  position  corresponding  to  the  maximum  critical  load  is 
ier  for  the  smaller  ratio  of  the  moduli.  For  E^/Et  =  0.1,  it  is  close 
1  ,r0  =  0.73.  Note  that  the  maxima  critical  loads  are  small  for  the 
ler  ratio  of  the  moduli,  too.  Furthermore,  the  bending  stiffness  of 
rir.g  has  a  stablizing  effect.  The  higher  the  stiffness,  the  higher  the 
ical  load.  Note  that  this  effect  is  more  pronounced  at  the  higher 
as  cf  rl/,ro  (r|/rQ  >  0.3).  On  figure  8,  one  can  also  see  curves  of 
:oai  lead  vs.  ring  location,  r^/rQ,  for  E2/Ei  =  1.0,  cir  =  0  and  various 
33  c-’  F~.  When  these  curves  are  compared  to  those  corresponding  to  olr 
',  they  depict  the  effect  of  neglecting  the  extensional  stiffness  of 
'■ing  (ore  HO]).  Even  for  this  low  value  of  aR  (=  0.10),  this  effect 
be  substantial  at  the  higher  values  of  both  r-\/rQ  and  Br.  Finally, 
tr.e  ring  is  located  at  the  outer  edge,  the  critical  load  does  not 
r.e  equal  to  4.2  D/r^  (for  E2/E^  =  1.0),  primarily  because  of  accounting 
r.e  axial  stiffness  of  the  ring  (see  equation  (7)). 

Eirure  i  depicts  the  same  results,  but  for  n  =  1  (asymmetric  mode), 

-  5.1  and  various  values  of  Yr  for  both  values  of  tne  ratio  of 
•*  ••  ’  =  1.0  and  0.1.  Clearly  the  loads  are  higher  than  these 
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"  ...  -.sing  o  -  J,  therefore,  buckling  is  governed  by  the  axisyrametric 
.  f  i  n  3 ly ,  it  .  j  shown  on  figure  9  that  the  higher  the  value  of  the 
t  r: ;  or.  el  stiffness  parameter ,  Y^ ,  the  higher  the  critical  load  (for  n  = 

s  xpected. 

11a  -  :~-^s  Multi-Annular  Plates 

n  „  s  r. r  example  consists  of  a  multi-annular  plate  clampeu  end  loaded 
■  t  •  uter  eage.  The  load  is  radial  and  axisymmetric .  The  parameters 
d-.-rc  :  m  this  case  are: 
v  =  0.3;  E2/E,  =  0,  0.1,  0.5,  and  1.0; 

f.r  n  ■  urj  =  0.03  and  =  1 .2;  and 

f:r  n  -  1,  as  =  0.0r'  and  6^  =  1.2  and  Yr  =  0.6 
T trends  which  are  observed  in  the  simply  supported  case  are  mostly  the 
;.:r  :.  The  results  are  presented  for  the  unstiffened  plate  on  figure 
ms  on  figures  11  and  12  for  the  stiffened  plate. 

I  comparison  between  the  unstiffened  case  and  the  stiffened  one  shows 
-  i  that  the  stiffener  increases  the  buckling  capacity  of  the  plate, 
r.  dices  of  the  presence  of  stiffening,  when  the  ratio  of  the  moduli  is 
tuckling  is  triggered  by  the  central  part  and  the  critical 
1  .  id  it  decreasing  as  the  common  boundary  approaches  the  outer  edge. 

h  ;te  tr.at  on  figure  11,  results  corresponding  to  the  annular  plate, 

-  6,  are  only  shown  for  rj/rQ  <  0.55.  It  has  been  shown  [11]  that 
:  :r  r'  >  9.55  higher  (n  >  0)  buckling  modes  govern  and  thus  the  critical 
.  n  i.rrespcnds  to  higher  values  of  n. 

.  ..  pr-.-ser.co  of  stiffening,  figures  11  and  12,  shows  critical  loads 
"  ‘  '  ir:r'-‘1  ir  P1  ites  ^ 11 2  /  E  i  =  0).  From  figure  12,  one  observes 

'-5  gov-.-rned  by  the  axisymmetric  mode.  This  is  true  for  the 
■ — j *  ■.  s  r  .  /  r  q  - ) .  ( . 
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1 
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Meal  loada  for  a  ring-st if fened  plate  (clamped;  n 


xirle  inner  plates,  E2/E'i  *  0,  but  with  0  <  E2/El  -  1-0,  and  a 
at  tr.e  common  joint,  buckling  is  governed  by  the  a  x  i  s  ymme  t  r  i  c 


■a  ^ 


D 


2 1  0 


■f  a  r.ing-stiffened  Plate  Using  Plate  Theory. 

us  section  a  different  approach  is  used  for  the  analysis  of  a 
us  circular  plate  stiffened  by  a  ring  somewhere  between  the  center 
ge .  The  stiffener,  in  this  case,  is  modelled  by  a  very  narrow 
late,  with  different  geometrical  and  material  properties.  The 
:arried  out  is  that  of  a  three-part  plate,  where  the  outer  and 
arts  have  the  same  properties,  and  the  middle  part  represents  the 


ner . 

nis  type  of  analysis  is  correct,  provided  that  plate  theory 
tior.s  hold  true  for  the  ring  part,  as  well.  The  ring  plate  geometry 
•operties  were  adjusted  in  order  to  yield  the  following  values  for 
ype  of  stiffnesses.  (This  was  done  in  order  to  be  able  to  compare 
esent  results  with  those  using  ring  stiffeners). 


a:.  =  0.1  with  £a  =  10 


0.2  with  =  20 


for  n  =  0 


a;_  ---  0.1  with  £3  =  1.0,  Yr  =  0.35  for  n  =  1. 


shown  graphically  on  figure  13  for  n  =  0  (axisymmetric )  and 


Q?  i 


for  n  =  1  . 


It  is  seen  from  these  figures  that,  the  comparison  is  fairly  good. 
1  ne  rust  realize  that  in  one  case  the  stiffener  is  considered  as  a  curved 
be  and  in  the  other  as  an  annular  plate. 

CONCLUSIONS 

Among  the  many  conclusions  of  the  parametric  studies  one  may  list  the 
f... owing  important  one s: 

Stiffening,  regardless  of  boundary  conditions  at  the  outer  loaded 
.  h.;5  a  staoilizing  effect.  The  higher  the  bending  stiffness  of  the 

the  signer  the  critical  load.  Moreover  the  extensional  ring 
bt::  fr  om:  improves  the  critical  load.  The  effect  of  bending  stiffness  is 
t:  for  ring  locations  higher  than  0.3  r0.  Similarly,  the  effect  of 

•  n  -.icnul  ring  stiffness  becomes  important  for  ring  locations  higher  than 
„ .  5  r ;  .  Tnese  conclusions  are  true  for  simply-supported  boundaries,  but 
one  trends  are  also  true  for  clamped  boundaries. 

i)  Modelling  the  stiffener  as  an  annular  plate  yields  critical  loads, 
w  ..-h  are  close  to  those  corresponding  to  modelling  the  ring  as  a  curved 
■  N^ts  that  the  latter  modelling  is  more  general,  because  it  can 

.  . ommodate  various  ring  geometries  instead  of  only  rectangular  cross 
a  e  :  a  i  t r  i  n  gs . 

i/  Trie  use  of  a  single  ring,  for  stiffening  a  circular  plate,  implies 
ilta  an  optimum  position.  For  a  homogeneous  plate,  the  optimum  position 
cm  men  ponds  to  approximately  r5/r0  =  0.7  for  simply  supported  boundaries 
r’  ro  -  0.5  for  clamped  ones.  These  optimal  locations  change  when  the 
aw  P ; r t s  of  the  plate  have  different  properties. 
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BUCKLING  OF  LAMINATED  CIRCULAR  PLATES 
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ABSTRACT 

Ihe  manuscript  deals  with  buckling  of 
s;..  et  '  anj  annular  circular  plates,  Bade  of 
■v  r  l  C  laminates  with  general  orientation  of 

iamlr.i  fibers.  The  plate  Is  subjected  to  an 
lr  ;.a-.e,  symmetric,  destabilising  load.  The  load 
is  Trilled  at  the  outer  edge  for  the  oomplete 
7  1  s *  e  and  a*,  the  outer  and  Inner  edges  for  the 
'•  rj.ar  one.  The  plate  Is  supported  In  various 
-T)S  along  its  t/sundarles. 

The  i ay e r  includes  the  derivation  of  the 
-over:.,  n  g  equations  and  discusses  the  various 
oft-. is  of  solution.  Numerical  results  are 
-noscmed  f  jr  a  clamped  plate. 

Introduct Ion 


F  1  *  r  reinforced  laminated  composites  are 
fairly  .’ini  widely  replacing  metals  as  construction 
tr  f i  lals,  especially  In  the  aanufactura  of 
surfs'-,  air  and  space  vehicles.  Thin  plates  and 
shells  are  popular  structural  elements.  Such 
nM  gur  at  lens  ,  in  service,  are  mostly  and  often 
subject*  l  to  external  loads  which  Bay  Induce 
:.  u  •'  k  ;  l  n  g  .  In  many  situations  buokllng  Is  an 
■u  r.  de  S  .  r  a  D  1  e  phenoaenon  and  It  forms  a  prlaary 
cons. deration  for  design.  Because  of  this, 
stability  of  laminated  plates  and  shells  continues 
to  attract  the  attention  of  the  structural 
engineer*  . 

M-st  of  the  analyses  reported  In  the  opsn 
literature  [  I  -8  ]  for  circular  configurations  are 
limited  to  polar  orthotropy,  a  case  Tor  which  th« 
varlo-.d  rigidities  are  constant  with  respect  to 
the  radial  and  c 1 rcunferent 1  el  coordinates. 

Another  group  of  plates  eaploya  rectilinear 
orthotropy  for  which  results  are  reported  In  tha 
open  literature  [6-7].  The  analysis,  for  this 
group  is  limited  to  vibration  probleaa  only.  Tha 
b  j  r  k  1 1  r  g  of  laminated  composite  plates  la  usually 
limited  to  rectangular  plates  [8-9].  A  very  de¬ 
tailed  survey  based  on  hundreds  of  references  by 
.el  [10],  on  composite  plates,  discusses 
results  and  analyses  for  only  rectangular 
C-rif  if  ••  all  ons  . 

The  manuscript  Includes  the  derivation  of  the 
governing  equations  and  discusses  the  various 
n-tnod')  of  solutions.  The  analysis  deals  with 
bu'-vllng  of  circular  laminated  ooaposlte  platea, 
loaded  by  ln-plane  aKlsymastrlo  forcea  at  the 
edges.  The  geometry  Is  limited  to  symmetric 
regular  angle  ply  configurations  and  lamination 
theory  [8,11]  Is  employed.  The  solution  proposed 
In  the  paper  Is  based  on  the  Galerkln  ana  the 
modified  Galerkln  procedures.  The  plates 
•o’,  aldered  herein  are  either  clamped  or 
s  1  mpl  y  sut  ;>or ted  . 
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Mathematical  Formulation 


The  geometry  consists  of  a  circular  laminated 
plate  subjected  to  axlsymmetrtc  loading  at  the 
edges  (see  Fig.  1).  The  laminate  Is  a  symmetric, 
regular,  angle-ply  one  [8]. 

The  governing  equations  for  the  primary  state 
are  [II]: 
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Moreover,  the  buckling  equation.  In  terms  of 
moments  and  foroe  resultants,  after  neglecting  the 
small  non-llnear  terms.  Is  (see  Appendix  B). 
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where  N°r>  N°#,  N°9  are  the  prlaary  state  stress 

resultants ,  M  M  H  a  are  the  bending  and 
rr  0  8  r  0 


twisting  Boaent  resultants;  w  the  transverse  dis¬ 
placement,  and  r,  0  the  coordinates  in  the  radial 
and  c  1  rcuaferent  lal  directions,  respectively. 

The  ln-plane,  primary  state,  kinematic 
relations  are  given  by 
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where  u,v  are  the  deformations  In  radial  and 
clrcumf er ent l al  directions,  respectively.  Note 
that  the  transverse  displacement  (w°)  la  zero  In 
the  primary  state.  The  relations  between  the 
stress  resultants  and  the  ln-plane  (reference 
plane)  strains  are 
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Substitution  of  Eqs.  (7)  Into  Eq.  (8)  ana  use  of 
Eqs .  (1)  yield 
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defined  1 ;;  Appendix  A. 

T n*»  tero3  Ajj  are  0-dependent  and  they  also 
Ceren.-j  on  the  relative  fiber  orientation  between 
tf  e  variola  layers . 

Two  formulations  can  be  used  for  the  primary 
stale  aial/sls.  Note  that  In  order  to  solve  the 
bur*  log  e  ju.it  l  >n  we  nee  1  the  pr  lnary  state  stress 


resultants. 
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applied  loadings.  If  some  of  the  l n- plana 
'■  -  i'Ujr  )  conditions  are  kinematic  then  the 
vloenatl  :  apji-.a-h  la  used  for  solving  the  primary 
*>  t  a  t  e  .  r.  the  other  hand.  If  the  l  n  -  pi  an® 
ary  -  *"11  Ilona  are  natural  then  the  foroe 
-r.  is  *)yod.  For  details  see  Appendix  B. 

'<*  -Inenjtl  ,  formulation  in  terns  of  u  and  v 
j’  Jthtf  upward.  The  governing  equations  are 
1  :/  substitution  of  k.  qs .  (3)  and  (4)  into 
.  /  1  e  i  -n  r.  g  : 
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where  (a^]  -  [Atj]  1  and  Lp  Is  the  differential 
operator  of  the  equation. 

The  coefficients  In  the  above  equation  are  In 
general  0-dependent.  There  Is  no  closed  form 
solution  and  an  approximate  one  can  be  derived 
through  the  Calerkln  procedure  by  assuming 


F(r , 0)  -  Fq ( r , 9 )  ♦  l  anFn(r,e)  (10) 

n*  1 

where  F  (r,9)  Is  chosen  in  such  a  way  that  It 
satisfies  the  boundary  conditions.  Thus,  the 
stresses  due  to  function  Fn(r,e)  must  vanish  at 
the  boundaries.  The  unknown  coefficients,  an,  are 
evaluated  through  the  Calerkln  Integrals,  or 
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derivatives  with  respect  to  0  are  not  null  and 
ro  ;  3  t  be  Included  In  the  equations. 

The  force  approach  determines  the  stress 
sir  l  it  1  on  by  the  Introduction  of  the  Airy 
stress  function  through  the  ln-plane  compatibility 
equation.  The  Airy  stress  function  and  primary 
stress  resultants  are  given  by  [12}. 
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and  the  compatibility  equation  Is: 
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The  functions  Fn  (r,9)  must  satisfy  the  following 
boundary  conditions; 
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at  r  -  R 


If  F  n ( r  ,  8  )  -  rn(r  )  cos  n«  then  the  above 

oondltlons  can  be  replaced  byi 
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Namely,  the  functions  fn(r)  can  be  chosen  to  be 
the  deflection  functions  of  a  clamped  plate. 

In  the  case  of  a  circular  plate,  loaded  at 
the  outer  edge  the  functions  arei 
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where  :•  -  Is  the  outer  radius.  , 

In  ease  of  an  annular  plate,  loaded  et 

'  r»  ,"ier  and  at  the  outer  edges,  tha  functions 

n s  .  :crcJ  are: 
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where  f*0  and  R|  are  the  outer  and  Inner  radii, 

r e a  r*  i  y  . 

r  n  tre  case  of  a  circular  plate,  loaded 
a  1 1  ■  i  r.tielr  i  ca  1 1  y  at  the  edge,  or  an  annular  plate 
loaded  at  Ida  outer  and  Inner  edges  with  the  saae 
Intensity,  an  analytical  closed  foro  solution 
exists,  provided  that  tho  plate  has  cylindrical 
o-ttiir:;,.  Tno  function  Is  given  by  F  -  Cr2, 
u r.ery  C  is  delcralned  througn  the  force  boundary 
coniltld's.  In  the  caee  of  a  laalnated  plate, 
alto  ;-?)-/• 30° >n3,  this  solution  la  aoceptable 
witnl-i  arse  er.glnerlng  accuracy.  The  error  In  the 
.. .'  c  pa  1 1  h  1  1 1 1  y  equation,  Eq.  (9)  Is  of  the  order 
12'’  ,  instead  cf  zero.  In  general,  the  solution 
;rz' delure  proposed  earlier  should  be  followed. 

Th  ;  buckling  equation,  in  tersa  of  the 
trjrsverse  displacement.  Is  derived  by  using  the 
roil '...eg  wlneaatlc  relations  [11): 
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W .  ,  i  .  i  are  the  radial,  circumferential 
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and  torsional  curvature  changes  of  the  plate 
mldsurfise.  The  aoa e n t - s t r a  1 n  relations  ere 
derived  tr, rough  Integration  through  the  thickness 
and  are  given  by  (Rote  that  there  Is  no  coupling 
between  bending  and  stretching) 
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:  : tut  ten  of  Eqs .  (IS)  and  (16)  Into  Eq.  (2) 
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where  N°  ,  N°  end  N°  are  the  radial,  olrcumfer- 
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entlal  and  shear  primary  stress  resultants.  The 
coefficients  In  the  above  equation,  Eq.  (16),  are 
8-dependent  and  the  partial  derivatives  with 
respect  to  0  must  bo  considered. 

Tho  solution  In  general,  even  for  the 
simplest  or  cases,  Is  very  complicated.  An 
approximate  solution  approach,  based  on  the 
Calerkln  procedure.  Is  employed.  It  Is  outlined, 
wl  .h  sufficient  detail.  In  the  next  section. 

Solution  of  the  Buckling  Equation 


Two  variations  of  the  Calerkln  procedure  are 
employed,  herein,  In  order  to  achieve  a  solution 
for  the  buckling  equation. 

The  first  one  Is  the  strtot  Calerkln 
procedure,  which  employs  s  series  solution 
(approximation)  such  that  each  term  In  the  series 
satisfies  all  boundary  conditions,  regardless  of 
their  nature  (kinematic  and/or  natural).  The 
second  one  la  the  modified  Calerkln  procedure  [13J. 
In  this  oase,  eaoh  term  of  the  series  need  only 
satisfy  the  klnemstlo  boundary  conditions.  Then 
minimization  of  the  error  Includes  boundary  terms 
along  with  the  Calerkln  Integrals. 

The  approximate  solution  Is  equal  to 
(trunoated  series): 
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where:  w  (r)  are  functions  of  r ;  A  undetermined 
oonatant  2nd  n  denotes  the  number  St  terms.  The 
assumed  function  to  be  used  with  the  regular 
Calerkln  procedure  [13),  must  satisfy  all  boundary 
conditions  et  the  various  edges  of  the  plate.  The 
radial  funotion  coefficients  ere  determined  In 
such  a  way  that  all  ths  boundary  conditions  are 
satisfied. 


I  to  clarify  this,  let  us  consider  the 

f  "  i  1  ov*  1  two  cases  for  a  coasclete  circular  plate. 

^  A  -la.nped  plate  at  r  -  R0, 

r.'\e  o.iodary  conditions  for  this, case  are: 
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Thus,  s ut  ) \ ; t ut I  on  cf  the  assumed  solution,  Eq . 
t  1  9  •  lr/.-i  sp  houna\ry  conditions  yields: 
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*n  l  *  !•*  t  r  I  ?  e  .1  through  the  application  of  the 
•'jale'-mn  prucelure  on  the  tuck  ling  equation,  Eq. 
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Tv-  >  conditions  for  this  plate  are  given 
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wr,er»  sn  is  determined  from  the  condition  that  the 
determinant  cf  the  coefficients  must  vanish  and  So 
t:  a  j  arc  considered  to  be  the  elements  of  an 
eigenvector.  Use  of  the  Galerkln  procedure  leads 
to  the  vanishing  of  n  Galerkln  Integrals,  or; 
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N  te  tna\  use  of  the  Galerkln  procedure  with  an 
ass  imed  solution  which  satisfies  all  the  boundary 
conditions  Is  advantageous  only  in  the  case  where 
all  boundary  conditions  are  kinematic.  In  case  of 
natural  toundary  conditions,  the  assumed  functions 
n  in  te  much  more  complicated  and  the  procedure 
m;g-it  ’■*  very  te.llOUS. 

A n  alternate  approach  Is  offered,  which  Is 
based  on  the  modified  Galerkln  procedure  as 
appears  1 a  [13].  In  this  case,  the  assumed 
solution  must  satisfy  only  the  kinematic  boundary 
r  editions,  and  the  unknowns  are  determined 
tnrough  the  following  modified  Galerkln  Integrals! 
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M  ,  g  n  the  raJlal  bending  eoaent  and  shear 
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*  A 

f  *  *•  i*»»  .^r  ■  i  red  with  the  approilMto  solutions! 

*  .  *  ,ir  e  t  f  e  r  ad  l  a l  bending  eomont  and  Shear 

f  a  n  8,  ;  ,  i*l  »i  the  fjoun  larles  denoted  by  the  l. 

N,j»er  leal  ('psuUj 

*  •  . liter  .ode  has  been  developed  for 
')  Mil  res-jlts.  The  code  eaploys 
’ .  .  ’■  .  e  ;  '  1}  ,  i  s  l  in  for  o<>th  the  [T  liar  y  state 
t  »./■*.  9  .  w»i!  as  for  the  b«j'**llng  analysis. 
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extension.  The  plate  Is  made  of  2m  plies,  with 
fiber  orientations  (0°/*  l80°/m)9,  where  m  is  the 
number  of  plies.  Note  that  when  a  -  3  (three 
piles),  the  configuration  Is  (  0°/ 60°/ - 60°  )  3 ;  when 
m  -  1,  the  configuration  Is  (0°/ k507"k5°/ 15° ) ;  and 
when  m  -  6,  it  finally  becomes  ( 0°/ 30°/ - 30°/ 30°/ 
-30°/30°)8.  All  piles  are  made  of  the  same 
material.  Graphl 1 tefEpoxy  (T300-5208),  with  7GJ 
fiber  volume  and  the  following  properties 
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The  buckling  load 

Is  nondlmenslonal 1  zed  by 

3  2 

tha  quantity  E^h  /I2  Rq(  1  -  v)2v2|h  wblch  Is 
common  for  all  configurations.  Thus, 

N  -  k  E,,h3/  12  R2  (1  -  v.-v... )  (29) 

rr  cr  1 1  o  12  21 

cr 

Note  that  h  and  R0  are  the  plate  thickness 
and  radius,  respectl  vely. 

The  primary  state  Is  solved  for  by  using 
the  procedure,  outlined  previously,  and  the 
functions  appearing  In  Eqs.  (13)  and  (Ik).  The 
Galerkln  integrals  are  evaluated  by  using 
numerical  Integration  with  Romberg’s  extrapol  at  I  on 
method.  The  allowable  error  Is  smaller  than  10"^. 

The  approximating  aeries  Involves  functions 
which  sre  the  buckling  modes  of  the  corre»|>ondlng 
Isotropic  problem. 

N 

W  (r,6)  -  y  (J  (a  r)  •  C  rn)  cos  n  e 
a  *•  n  n  2n 

n-0 

where,  Jn(  )  are  Bessel  functions  of  the  first 
kind  of  order  m  an  13  th8  square  root  of  the 
buckling  coefficient  or  the  lsotrfplc  critical 
stress  resultant  (corresponding  to  the  nth  mode; 

If  n  •  0  an  •  Ik. 68):  and  C2n  ar*  undetermined 
constants,  evaluated  by  satisfaction  of  the 
boundary  conditions. 

The  limited  results  are  shown  graphically 
on  Figs.  2  and  3.  Fig.  2  shows  a  plot  of  kcr 
versus  m,  which  Implies  various  quas l -  1  sotropi c 
(In  extension)  constructions.  It  Is  seen  from 
this  figure,  that  kcr  Is  independent  of  m.  This 
Is  so  because  the  average  flexural  stiffnesses, 
primarily  D)j  &  D2z>  are  the  same  for  all  m. 
Next,  the  configuration  corresponding  to  m  ■  3  was 
taken  and  E??/Eti  “a»  arbitrarily  changed  rroro  the 
correot  value  of  0.03665  to  the  value  of  one 
(almost  Isotropic)  by  keeping  the  other  values 
constant.  It  is  seen  from  Fig.  3  that  the  value 
of  kcr  varies  from  6.(6  to  12.3'.  Note  that  the 
value  of  kC(.  yor  pne  Isotropic  construction  Is 
Ik. 68.  The  difference  between  the  two  values  Is 
attributed  to  the  effect  of  the  shear  stiffness. 
On  Fig.  3,  and  for  E22/E||  -  1,  two  additional 
values  oT  kCr  ar8  9hown,  In  order  to  depict  the 
G-sfrect . 
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Append  1  «  A  :  Compliance  Matrix  Coefficients  In  6 

Coordinates 

This  appendix  outlines  the  procedure  for 
determining  the  coefficients  of  the  aatrlx  In 
terms  of  the  fiber  orientation  angle  of  eeoh 
laminate  relative  to  the  0-0  ails. 
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Rotation  of  th«  coapllance  matrix,  with  respect  to 
van  angle  (a~e)f  yields  the  following  coefficients: 
wh«re» 
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X  v 

where  o  -  cos  (a,  -  0)  end  a  -  Ki  (o  -  0)  \ 

k  X  k  N 

X  x 

Thus,  In  the  r,0-ooor dl net e  ayatfcis  the  stress'' 
atete  Is  related  to  the  state  of  strain  by.  X 

X.  \ 
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Append  1 x  Bt  Derivation  of  the  Governing  Equations 


This  appendix  la  devoted  to  a  detailed 
derivation  of  the  governing  equations.  The 
positive  directions  of  the  applied  end  Internal 
forces  ars  desorlbed  In  Fig.  B.l. 


The  stress-atrsln  relation  for  each  lamina  wider 

consideration  la: 
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(a)  In-plane  Resultants  (b)  Internal  Bending 

Forces  end  Torolon  Momenta 
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(c>  Internal  Vertical  (d)  Positive  Deforest  Ion 
Shear  Forces 


Fig.  B.l  Internal  Forces  A  Deforaatlon  of 
Differential  Area  of  the  Plate 
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The  In-plane  oqulllbrlua  equations  In  l eras 
of  the  Internal  forces  are: 
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and  tne  transveree  are  In  teras  of  lnplane  foroes 
and  Internal  rooaents,  using  the  equations  above. 
Is 
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The  bending  aoaent  end  the  vertical  deflection  In 
the  pr leery  state  designated  with  (  )°  are  aero. 

The  additional  straaa,  la  ohoaan  to  be  very 

saall  coapared  to  N°ji  thus,  after  substituting  It 
In  Eqs.  (B.l)  and  (8.2)  yields: 
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and  the  buckling  equation  In  the  vertical 
direction  Is: 
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Slnoe  Hjj  <  <  11^  ,  neglect  of  these  teras  and 

oaltttng  the  •  sign  of  the  deflection  and  the 
aoaenta  yields: 
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Tne  curvatures  are: 
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Using  ns  ccimon  assiaaptlon  In  plats  theory  yields 
tr.e  following  relation  between  the  Internal  forces 
and  tne  strains  and  curvatures: 
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Since  the  various  stiffnesses  are  8-dependent,  the 
derivative  of  the  coefficients  with  respect  to  e 
are  different  fro«  aero  and  Bust  be  Included  In 
the  equations. 


The  In-plane  stress  distribution  can  be 
derived  using  the  Airy  function,  through  the 
coapablllty  equation.  The  Airy  function  and 
stress  relations  are: 
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The  foreulatlon  using  deformations  u,v  and  w 
yields  the  following  governing  equations: 


The  ooapatlblllty  equation  Is: 
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Substltuton  In  Eq.  (8. 19)  and  use  of  Eqs.  (8.6) 
and  (B. 1 3>  yields: 
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«r,*re  Lp  is  tne  differential  operator  described  In 

c fi*  e  4ua t ion . 

Tr.e  buckling  equation  is  given  by: 
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Fig.  1  toad,  boundary  conditions  and  geoaetry  of 
a  typical  laalnatad  circular  plate 


